A re-activation model for 169Yb intensity modulated brachytherapy sources accounting for spatiotemporal isotopic composition.

BACKGROUND: Intensity modulated brachytherapy based on partially shielded intracavitary and interstitial applicators is possible with a cost-effective 169Yb production method. 169Yb is a traditionally expensive isotope suitable for this purpose, with an average γ-ray energy of 93 keV. Re-activating a single 169Yb source multiple times in a nuclear reactor between clinical uses was shown to theoretically reduce cost by approximately 75% relative to conventional single-activation sources. With re-activation, substantial spatiotemporal variation in isotopic source composition is expected between activations via 168Yb burnup and 169Yb decay, resulting in time dependent neutron transmission, precursor usage, and reactor time needed per re-activation. PURPOSE: To introduce a generalized model of radioactive source production that accounts for spatiotemporal variation in isotopic source composition to improve the efficiency estimate of the 169Yb production process, with and without re-activation. METHODS AND MATERIALS: A time-dependent thermal neutron transport, isotope transmutation, and decay model was developed. Thermal neutron flux within partitioned sub-volumes of a cylindrical active source was calculated by raytracing through the spatiotemporal dependent isotopic composition throughout the source, accounting for thermal neutron attenuation along each ray. The model was benchmarked, generalized, and applied to a variety of active source dimensions with radii ranging from 0.4 to 1.0 mm, lengths from 2.5 to 10.5 mm, and volumes from 0.31 to 7.85 mm3, at thermal neutron fluxes from 1 × 1014 to 1 × 1015 n cm-2 s-1. The 168Yb-Yb2O3 density was 8.5 g cm-3 with 82% 168Yb-enrichment. As an example, a reference re-activatable 169Yb active source (RRS) constructed of 82%-enriched 168Yb-Yb2O3 precursor was modeled, with 0.6 mm diameter, 10.5 mm length, 3 mm3 volume, 8.5 g cm-3 density, and a thermal neutron activation flux of 4 × 1014 neutrons cm-2 s-1. RESULTS: The average clinical 169Yb activity for a 0.99 versus 0.31 mm3 source dropped from 20.1 to 7.5 Ci for a 4 × 1014 n cm-2 s-1 activation flux and from 20.9 to 8.7 Ci for a 1 × 1015 n cm-2 s-1 activation flux. For thermal neutron fluxes ≥2 × 1014 n cm-2 s-1, total precursor and reactor time per clinic-year were maximized at a source volume of 0.99 mm3 and reached a near minimum at 3 mm3. When the spatiotemporal isotopic composition effect was accounted for, average thermal neutron transmission increased over RRS lifetime from 23.6% to 55.9%. A 28% reduction (42.5 days to 30.6 days) in the reactor time needed per clinic-year for the RRS is predicted relative to a model that does not account for spatiotemporal isotopic composition effects. CONCLUSIONS: Accounting for spatiotemporal isotopic composition effects within the RRS results in a 28% reduction in the reactor time per clinic-year relative to the case in which such changes are not accounted for. Smaller volume sources had a disadvantage in that average clinical 169Yb activity decreased substantially below 20 Ci for source volumes under 1 mm3. Increasing source volume above 3 mm3 adds little value in precursor and reactor time savings and has a geometric disadvantage.

The population percentile allowance method for determining systematic spatial error tolerances for temporary intensity modulated brachytherapy.

BACKGROUND: Multiple approaches are under development for delivering temporary intensity modulated brachytherapy (IMBT) using partially shielded applicators wherein the delivered dose distributions are sensitive to spatial uncertainties in both the applicator position and shield orientation, rather than only applicator position as with conventional high-dose-rate brachytherapy (HDR-BT). Sensitivity analyses to spatial uncertainties have been reported as components of publications on these emerging technologies, however, a generalized framework for the rigorous determination of the spatial uncertainty tolerances of dose-volume parameters is needed. PURPOSE: To derive and present the population percentile allowance (PPA) method, a generalized mathematical and statistical framework to evaluate the tolerance of temporary IMBT approaches to spatial uncertainties in applicator position and shield orientation. METHODS: A mathematical formalism describing geometric applicator position and shield orientation shifts was derived that supports straight and curved applicators and applies to serial and helical rotating shield brachytherapy (RSBT) and direction modulated brachytherapy (DMBT). The PPA method entails defining the percentage of a patient population receiving a given therapy that is, allowed to receive dose-volume errors in the target volume and specified organs at risk of a defined percentage or less, then determining what combinations of applicator position and shield orientation systematic errors would be expected to produce that outcome in the population. The PPA method was applied to the use case of multi-shield helical 169 Yb-based RSBT for cervical cancer, with 45° and 180° shield emission angles. A total of 37 cervical cancer patients were considered in the population, with average (± 1 standard deviation) HR-CTV volumes of 79 cm3  ± 37 cm3 and optimized baseline treatment plans (no spatial uncertainties applied) created for each patient to meet dose-volume requirements of 85 GyEQD2 (equivalent uniform dose in 2 Gy fraction), with D2cc tolerance doses of 90 GyEQD2 , 75 GyEQD2 , and 75 GyEQD2 for bladder, rectum, and sigmoid colon, respectively. RESULTS: For the PPA requirement that 90% of cervical cancer patients receiving multi-shield helical RSBT could have a maximum dose-volume uncertainty of 10% for high-risk clinical target volume (HR-CTV) D90 (minimum dose to hottest 90%) and bladder, rectum, and sigmoid colon D2cc (minimum dose to hottest 2 cm3 ), the tolerance systematic applicator position and shield orientation uncertainties were approximately ± 1.0 mm and ± 4.25°, respectively. For ± 1.5 mm and ± 5° systematic applicator position and shield orientation tolerances, 90% of the patients considered would have a maximum dose-volume uncertainty of 12.8% or less. CONCLUSION: The PPA method was formalized to determine the temporary IMBT spatial uncertainty tolerances that would be expected to result in an allowed percentage of a population of patients receiving relative dose-volume errors above a defined percentage. Multi-shield, helical 169 Yb-based RSBT for cervical cancer was evaluated and tolerances determined, which, if applied on each treatment fraction, would represent an extreme situation. The PPA method is applicable to a variety of temporary IMBT approaches and can be used to rigorously determine the design parameters for the delivery systems such as mechanical driver motor accuracy, shield angle backlash, applicator rotation, and applicator fixation stability.

Efficient 169 Yb high-dose-rate brachytherapy source production using reactivation.

PURPOSE: To present and quantify the effectiveness of a method for the efficient production of 169 Yb high-dose-rate brachytherapy sources with 27 Ci activity upon clinical delivery, which have about the same dose rate in water at 1 cm from the source center as 10 Ci 192 Ir sources. MATERIALS: A theoretical framework for 169 Yb source activation and reactivation using thermal neutrons in a research reactor and 168 Yb-Yb2 O3 precursor is derived and benchmarked against published data. The model is dependent primarily on precursor 168 Yb enrichment percentage, active source volume of the active element, and average thermal neutron flux within the active source. RESULTS: Efficiency gains in 169 Yb source production are achievable through reactivation, and the gains increase with active source volume. For an average thermal neutron flux within the active source of 1 × 1014  n cm-2  s-1 , increasing the active source volume from 1 to 3 mm3 decreased reactor-days needed to generate one clinic-year of 169 Yb from 256 days yr-1 to 59 days yr-1 , and 82%-enriched precursor dropped from 80 mg yr-1 to 21 mg yr-1 . A resource reduction of 74%-77% is predicted for an active source volume increase from 1 to 3 mm3 . CONCLUSIONS: Dramatic cost savings are achievable in 169 Yb source production costs through reactivation if active sources larger than 1 mm3 are used.

Interstitial rotating shield brachytherapy for prostate cancer.

PURPOSE: To present a novel needle, catheter, and radiation source system for interstitial rotating shield brachytherapy (I-RSBT) of the prostate. I-RSBT is a promising technique for reducing urethra, rectum, and bladder dose relative to conventional interstitial high-dose-rate brachytherapy (HDR-BT). METHODS: A wire-mounted 62 GBq(153)Gd source is proposed with an encapsulated diameter of 0.59 mm, active diameter of 0.44 mm, and active length of 10 mm. A concept model I-RSBT needle/catheter pair was constructed using concentric 50 and 75 µm thick nickel-titanium alloy (nitinol) tubes. The needle is 16-gauge (1.651 mm) in outer diameter and the catheter contains a 535 µm thick platinum shield. I-RSBT and conventional HDR-BT treatment plans for a prostate cancer patient were generated based on Monte Carlo dose calculations. In order to minimize urethral dose, urethral dose gradient volumes within 0-5 mm of the urethra surface were allowed to receive doses less than the prescribed dose of 100%. RESULTS: The platinum shield reduced the dose rate on the shielded side of the source at 1 cm off-axis to 6.4% of the dose rate on the unshielded side. For the case considered, for the same minimum dose to the hottest 98% of the clinical target volume (D(98%)), I-RSBT reduced urethral D(0.1cc) below that of conventional HDR-BT by 29%, 33%, 38%, and 44% for urethral dose gradient volumes within 0, 1, 3, and 5 mm of the urethra surface, respectively. Percentages are expressed relative to the prescription dose of 100%. For the case considered, for the same urethral dose gradient volumes, rectum D(1cc) was reduced by 7%, 6%, 6%, and 6%, respectively, and bladder D(1cc) was reduced by 4%, 5%, 5%, and 6%, respectively. Treatment time to deliver 20 Gy with I-RSBT was 154 min with ten 62 GBq (153)Gd sources. CONCLUSIONS: For the case considered, the proposed(153)Gd-based I-RSBT system has the potential to lower the urethral dose relative to HDR-BT by 29%-44% if the clinician allows a urethral dose gradient volume of 0-5 mm around the urethra to receive a dose below the prescription. A multisource approach is necessary in order to deliver the proposed (153)Gd-based I-RSBT technique in reasonable treatment times.